Ring Exchange and Correlated Fermions 

Michel Roger 

DSM-DRECAM, Service de Physique de VEtat Condense 
Orme des Merisiers, CEA Saclay 
91191 Gif sur Yvette Cedex. France 



Abstract 

The concept of exchange in strongly-correlated fermions is reviewed with emphasis on the generalization of the 
Heisenberg pair exchange to higher order n-particle permutations. The "frustration" resulting from competing 
ferromagnetic three-spin exchange and antiferromagnetic two- and four-spin exchanges is illustrated on a two- 
dimensional model system: solid 3 He films. Recent experimental results proving the presence of significant four-spin 
exchange interactions in the CuCh plaquettes of high T c cuprates are reported. 



Key words: Quantum spin systems, 3 He, Cuprates. 



1. Historical introduction 

The concept of exchange interactions in almost- 
localized correlated fermions first appeared in the 
pioneering papers of Heisenberg [1] and was for- 
mulated in a more general mathematical way by 
Dirac[2,3]. Although the early papers by Dirac 
already contain the general expression of the 
Hamiltonian in terms of n-particle permutations, 
no attention was paid, up to the sixties, to higher 
orders than the "pair-exchange" Heisenberg term. 
Thouless[4] was the first to point out that higher 
order exchanges as three- and four-spin cyclic per- 
mutations might be important in quantum solids 
like 3 He. But only ten years later, striking exper- 
imental results on nuclear magnetism in the bcc 
phase of solid 3 He, in the millikelvin range[5,6,7], 
were interpreted by Hetherington, Delrieu and 
Roger through a ring-exchange model with two- 
three- and four-particle interactions of comparable 



magnitude[8,9,10]. From general considerations 
put forward by Thouless, permutations of even 
parity (like three-particle cycles) induce ferromag- 
netism while odd permutations (pair and four- 
particle exchanges) favor antiferromagnetism, 
and the striking phase diagram of bcc solid 3 He 
corresponds to a highly frustrated quantum-spin 
system with competing three- and four-spin per- 
mutations. During the last two decades a lot of 
progress has been accomplished in the investiga- 
tion of solid 3 He films adsorbed on graphite, a sim- 
ple model system exhibiting even more frustration 
since the frustrated nature of ring-exchange in- 
teractions is enhanced by the frustrated geometry 
of the triangular lattice. The conceptual beauty 
of solid 3 He lies in the fact that the system is 
simple enough (the pair interactions between 3 He 
atoms are mainly of hard core nature) to allow 
the calculation of exchange frequencies from first 
principles [11] and a quantitative comparison with 
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experimental results. 

The relevance of the multi-spin exchange con- 
cept is not restricted to the physics of nuclear 3 He 
spins. Delrieu[12] suggested that three-spin ex- 
change might be dominant in the two-dimensional 
electron Wigner solid near the quasi-classical 
limit, and this has been corroborated through 
WKB calculations[13]. More recent Monte-Carlo 
simulations have proved that competing three and 
four-spin exchange interaction should occur near 
melting [14]. 

The first fourth-order t/U expansion of the 
one-band Hubbard Hamiltonian in terms of four- 
spin interactions has been published in 1977 by 
Takahashi[15]. Soon after the discovery of high-T c 
superconductors, Roger and Delrieu[16] suggested, 
on the basis of an expansion of the three-band 
Hubbard model, that four-spin exchange might 
be significant in the CuC>2 planes of cuprates. 
During the last ten years many experimental 
results have revealed the presence of four-spin 
exchange interactions in the Cu-0 plaquettes of 
cuprates[17,18,19] and copper-based spin-ladder 
materials[20,21,22,23]. 



2. Dirac formalism with illustration on the 
one-band Hubbard model. 

Dirac formalism is introduced in the framework 
of degenerate perturbation theory. The Hamilto- 
nian is written H = Hq + V, where Hq describes 
independent particles and V is a perturbation. The 
"unperturbed" degenerate ground-states for N dis- 
tinguishable particles can be written as products 
of independent particle states: 



m = |ai} (1) |a2) (2) |a3> (3) ---|^) W 



(1) 



which means that the particle numbered (i) is in a 
state \cti). Each of these states is itself a product 
of two kets corresponding to the orbital and spin 
variables respectively: \a.i) = \Ri)\<Ji) 

^■■■\R n ) (N) \a N ) (N) 
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Hence \ip) = \4> R )\^' 7 ) appears as a product of an 
orbital wavefunction: 
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and a spin wavefunction: 
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Note that for the half- filled Hubbard model, \Ri) 
simply represents the site occupied by the parti- 
cle (i). The permutation P of two particles can be 
expressed as a product of two operators: 



p _ pHp<7 



(5) 



P R acting on orbitals and P a acting on spin vari- 
ables. If the Hamiltonian does not depend explic- 
itly on the spin, we can as a first step concentrate 
on the orbital part of the wave function and solve 
the eigenvalue problem: 



H\^ R )=E\4> R ) 



(6) 



for the orbital wavefunction \4> R ) describing N dis- 
tinguishable particles. The ground-state of the un- 
perturbed part Hq of the Hamiltonian is N! fold 
degenerated and the corresponding subspace Hq is 
spanned by the N! states: 



> = |^> (1) l^> (2) 



(7) 



where {yi, v^, • • ■ , vn} represents a permutation 
P of the N integers {1, 2, • • • , N}. We now apply 
degenerate perturbation theory[25,26] to the per- 
turbed Hamiltonian H = Hq + V. At first order, 
the splitting of the N! degenerated energy levels is 
given by the eigenvalues of the Hamiltonian 
defined by its matrix elements: 



v a,b — 



P R VP R \<P R ) 



(8) 



where P R and P R are two permutations of the 
symmetric group S n ■ Since V is invariant with re- 
spect to any permutation, we can write: 

v$ = (4> R \VP R P b R \4> R ) = (^ r \vp r \4> r ) = 

where P = P R P R 7 and the eigenvalue problem, 
restricted to the subspace i7o can be formally rep- 
resented by the Hamiltonian: 



h^ = - £ 



V { ^P R 



(9) 



where the summation runs over permutations P R 
of the symmetric group Sn- This result extends 
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straightforwardly to higher order degenerate per- 
turbation theory: the higher orders are expressed 
in terms powers of V and projection operator Pq 
on ft [25,26], and these operators commute with 
permutation operators. Hence, at arbitrary order 
in degenerate perturbation theory, we can write: 

Hks- V P pR ( 10 ) 

p R es N 

We now have to introduce the spin degrees of free- 
dom an express that the global wavefunction is 
completely antisymmetric. As a general result from 
Group theory a completely antisymmetric wave 
function can be expressed by the following bilinear 
expression [2 7]: 

iv>> = E c ^a>i^> (ii) 

where \4> R ) represents a linear combination of dif- 
ferent permutations P R \<fi R ) corresponding to a 
given irreducible representation of the symmetry 
group schematized by a Young diagram, while 
represents a linear combination of permuta- 
tion Pg corresponding to the representation 
associated with the "complementary" Young di- 
agram, obtained by exchanging the lines and the 
columns. For spin-1/2, a complete antisymmetri- 
sation of the spins cannot be realised over more 
than 2 variables, hence the corresponding Young 
diagrams have at most two lines and each diagram 
corresponds to a given value of the total spin 5*. 
It is then possible to establish a correspondence 
between the expression [Eq. (10)] of the Hamilto- 
nian acting only on the orbital variables with an 
equivalent Hamiltonian acting only on the spin 
variables. Expressing the antisymmetry of the 
wave function: 

P|V>> = (-l) p |V) (12) 

where p is the parity of the permutation, we can 
write: 

P R P°\i,) = (-lflV) (13) 
and multiplying to the left by (P^) -1 : 

p^) = {-iy{p R )-^) (w) 



Taking into account that in Eq. (10) P R and the 
inverse permutation (P^) -1 appear with the same 
weight Vp, the Hamiltonian is written equivalently 
in spin space: 

J] {-lfVpP a (15) 

P°es N 

where the sum is over all permutations of the sym- 
metric group Sm acting on spin variables. 

Large U expansions of the one-band Hubbard 
model[15] have generally been performed within 
the framework of second quantization (i.e working 
on the total antisymmetric wavefunction). Dirac 
formalism allows a simpler and more physical 
derivation. We start with N distinguishable parti- 
cles on N discrete sites, each site containing zero, 
one or two particles. There is an onsite Coulomb 
energy U for putting two particles on the same site, 
and each particle has a probability t to hop from 
one site to the nearest neighbor. The Hamiltonian 
does not depend explicitly on the spin variables. 

In perturbation theory, we start from the infi- 
nite U limit: each site is occupied by one particle. 
Since the particles are considered as distinguish- 
able, the ground state is N! fold degenerated. We 
apply degenerate perturbation theory to express 
the energy splitting as an expansion in t/U when 
U is large but finite[25,26]. Up to fourth order, the 
effective Hamiltonian restricted to the subspace fio 
spanned by the N! fold degenerated ground states 
is written[15, 25,26]: 

P VSVP +[PoVSVSVSVP -P VS 2 VP VSVP } 

where Pq is the projection operator onto the sub- 
space n and S = (1 - P )/(Eq - H ). This ex- 
pression can be expanded in term of two, three and 
four-particle permutations P R acting on coordi- 
nates: 

(1) (2) (3) 

h = -j 1 ^p^-j 2 Yp r -j 3 Yp r - 

<i,j> <i,k> <i,l> 

<i,j,k> 

k £ [p^i + ip^r 1 ] 

<i,j,k,l> 
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Fig. 1. Evaluation of the four-particle cyclic permutation 
frequency. Only the term P VSVSVSVP in the fourth 
order expansion contribute to the cyclic exchange of the 
four particle coordinates i —* j, j k, k —* I and ( — > i. 
Each of these four particle hops is accomplished by one of 
the four V in the previous term and at each step a factor 
—t a R I {vU) where v is the number of doubly occupied sites, 
in the corresponding intermediate virtual state, has to be 
taken into account (this factor arises from the operator S). 
There are 4!=24 terms depending on the order in which 
the four hops 1, 2, 3, 4, are successively realised. Eight 
of them [corresponding to the sequences (1324), (1342), 
(3124), (3142), (2413), (4213), (2431), (4231)] contribute 
to a factor — t 4 /(2U 3 ), the sixteen others contribute to 
-t 4 /U 3 . The total contribution is then -20t 4 /U 3 . 

with: Ji = 2(< 2 /[/)(l + At 2 /U), J 2 = 12i 4 /[/ 3 , 
J 3 = 2t 4 /U 3 , J T = lO^/U 3 and K = 20t 4 /C/ 3 for 
the square lattice. The sums in the first line are on 
first, second and third neighbor pairs respectively. 
The sum run over triangles in the second line and 
over square plaquettes in the last line. As a typical 
example, the Fig. 1 illustrates the evaluation of the 
four-particle cyclic permutation term. 

Since cyclic two and four-particle permutations 
are of odd parity, while three-particle permutations 
are even, we obtain the corresponding exchange 
Hamiltonian, with permutation operators P a act- 
ing on spin variables:: 

(1) (2) (3) 

Hex = j, J2 p?o + J* E p ik + Js E p u- 

<i,j> <i,k> <i,l> 

j t E i p ? j k + (p? jk )- 1 ] + 

<i,j,k> 

K E l P W + ( p ijkl)- 1 } (16) 

<i ,j,k,l> 

Using the identity: P?j = (1 + cr, • <Jj)/2 and the 
decomposition of any permutation in a product of 
transpositions, it is easy to check that this Eq. 16 is 
identical to the result given by Takahashi. However 



our previous expression clearly separate the ferro- 
magnetic contributions arising from three-particle 
exchange and the antiferromagnetic ones coming 
from pair- and four-particle exchange. The formu- 
lation in terms of permutation operators is also 
generally more convenient for the calculation of the 
thermodynamic properties of a n-particle exchange 
Hamiltonian. 



3. Cyclic exchange in solid 3 He and in the 
Wigner Solid from first principles 

Dirac formalism has been generalized by Thou- 
less for the exchange-problem in solid 3 He [4] and 
by Herring for many-electron systems [24]. This 
has been applied later to ab-initio calculations of 
multi-spin exchange frequencies in solid 3 He and 
in the two-dimensional Wigner-Solid. A pioneer- 
ing multidimensional WKB calculation of the ex- 
change frequencies has been applied to the two- 
dimensional Wigner solid near the quasi-classical 
limit [13]. It has proved the conjecture by Delrieu 
that three-particle exchange should dominate in 
that regime. Although, at physical densities, solid 
3 He is far from the quasi-classical regime, the same 
WKB approach was extrapolated with qualitative 
results concerning the hierarchy of various cyclic- 
exchange frequencies. 

With the rapid increase of the performances of 
computers, ab-initio calculations of ring-exchange 
frequencies through Path Integral Monte Carlo 
(PIMC) became possible. In contrast to WKB 
evaluation they are relevant to the whole range 
of physical densities, including solid phases near 
melting where quantum fluctuations are large. 
The exchange frequencies for bec solid 3 He[ll] 
near melting, are in excellent agreement with the 
experimental thermodynamic properties, which 
can be inferred through exact high-temperature 
series expansions. 

Although the long-range Coulomb potential is 
completely different from the hard core poten- 
tial between Helium atoms, WKB calculations 
near the classical limit[13,28] and PIMC in the 
whole density range[14] give striking analogies 
between solid 3 He films and the electron Wigner 
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solid, on the same triangular lattice. In both sys- 
tems, three-particle exchange dominates near the 
quasi-classical limit, leading to ferromagnetism, 
while near melting, antiferromagnetic two- and 
four-particle exchange compete with ferromag- 
netic three-particle exchange leading to highly 
frustrated antiferromagnetic systems. 

Close to melting, both systems might be pro- 
totypes of two-dimensional "spin liquids" [29] with 
no long range order at T=0. Experiments on the 
magnetism of the Wigner Solid are now at the very 
beginning [30], but the nuclear magnetism of solid 
3 He films has now been extensively studied. 



4. Two-dimensional solid 3 He: a model 
system with tunable frustration. 

Exfoliated graphite offers large (10x10 nm) flat 
crystallites with a strong adsorption potential for 
He. Up to two solid He films can be adsorbed on 
this substrate, within a wide density-range. The 
magnetism of 2D helium films has been thoroughly 
studied, down to 10-100 in: 

(i) submonolayer solid 3 He films on graphite[31] 

(ii) solid 3 He films on graphite coated by a high 
density 3 He[32,33,34,35,36] or 4 He[38,39] 
layer or by a HD bilayer[37]. 

The main difference between (i) and (ii) is that, due 
to a stronger adsorption potential, the first layer 
is closer to the pure 2D case. At the same density 
exchange frequencies in the first layer are 10 times 
lower than in the second layer where 3 He atom 
permutations are made easier through excursion 
in the third dimension. However the same general 
trends are observed: 

- A low-density commensurate solid phase is sta- 
bilized by the graphite potential in case (i) and 
by the periodic potential of the first layer in case 
(ii). Its susceptibility has an antiferromagnetic 
character with a negative Curie- Weiss constant. 
Both susceptibility and specific heat present 
anomalous deviations with respect to the usual 
asymptotic behavior at high temperature (Curie 
Weiss law for the susceptibility and T -2 behav- 
ior for the specific heat). These features are due 
to the frustrated nature of the system and can 



be quantitatively fit through the ring-exchange 
model[36]. Down to 10 /iK no long-range order- 
ing has been detected, which seems to confirm 
the presence of a spin-liquid state. An ultra-low 
temperature specific- heat measurement [35] and 
susceptibility measurements [38,39] put a higher 
limit of 10-100 /iK on a possible gap in the ex- 
citation spectrum. Hence solid 3 Hc films could 
be the first experimental realization of a gapless 
spin liquid. 

- At higher densities, an incommensurate solid 
phase is observed. The susceptibility has a pos- 
itive Curie Weiss constant, but the unusual 
features of the specific heat indicate that we 
have a "frustrated ferromagnetic" phase with 
relatively larger three-spin exchange but signifi- 
cant competing four-particle exchange [36]. The 
frustration is easily "tunable" in the sense that 
when the density is further increased, three-spin 
exchange dominates leading to a more conven- 
tional fcrromagnet. 
Some anomalies in the specific heat are present 
only in the commensurate phases. A very interest- 
ing interpretation has been recently proposed by 
H. Fukuyama[31], in terms of "ground-state vacan- 
cies" surrounded by magnetic polarons. This ex- 
citing idea established a close link with strongly 
interacting almost-localized electron systems close 
to a Mott-Hubbard transition (e.g the physics of 
cuprates!). 



5. Four-spin exchange in high-T c Cuprates 

In the insulating phases of Cuprates, the first 
experimental results in contradiction with a pure 
Heisenberg model were the anomalously broad 
multi-magnon Raman spectrum[17] and the in- 
frared optical absorption by phonon assisted 
multimagnon excitations [40] in La 2 Cu04. Both 
results can be understood with the occurrence of 
four-spin exchange [17, 20]. However the key exper- 
iments which really invalidates the pure Heisen- 
berg model is the inelastic neutron study of the 
magnon spectrum near the Brillouin zone [18]. It 
is in agreement with the presence of ferromagnetic 
three-spin exchange terms Jt as those appearing 
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in the expansion of the one-band Hubbard model 
(Eq. 16). However, due to accidental cancellations 
occurring for the Neel phase of the square lattice, 
the linear spin- wave spectrum is insensitive to the 
presence of four-spin exchange operators [19]. An 
excellent way to probe directly the presence of 
four-spin terms is the study of spin-spin correla- 
tions in the paramagnetic phase. These correlation 
have been recently measured through diffuse mag- 
netic scattering of polarized neutrons in La2Cu04 
at 350 K < T < 500 K[19}. The results are com- 
pared to predictions from high-temperature series 
expansions of the Hamiltonian. They prove defini- 
tively the presence of four-spin interaction with a 
ratio K/J\ of about 20%. 

The same Cu-0 plaquettes are present in copper 
based quasi-lD spin ladder materials Sri4Cu2404i 
and LagCa8Cu2404i, with approximately the same 
Cu-0 bonds and distances as in La2Cu04. The first 
neighbor pair exchanges along rungs and legs are 
expected to be roughly equal J rung ~ -heg — J- In 
that case the theory predict a gap A « 0.5 J in the 
excitation spectrum. The experimental value of the 
gap is substantially lower: A exp « 0.35J[21,22]. 
A relatively small four-spin exchange K/J « 20% 
accounts for the decrease of the gap [23]. This ef- 
fect is related to the frustration due to four-spin 
exchange, which increases the density of low lying 
eigenstates. 

The concept of multi-particle exchange, pio- 
neered by Dirac at the beginning of the 20th cen- 
tury, now appears as essential in the physics of 
strongly- correlated fermion systems which have 
raised considerable interest during the last decades. 
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